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ABSTRACT 



Context. The M-type star Gliese 581 is orbited by at least one terrestrial planet candidate in the habitable zone, i.e. GL 581 
d. Orbital simulations have shown that additional planets inside the habitable zone of GL 581 would be dynamically stable. 
Recently, two further planet candidates have been claimed, one of them in the habitable zone. 

Aims. In view of the ongoing search for planets around M stars which is expected to result in numerous detections of poten- 
tially habitable Super-Earths, we take the GL 581 system as an example to investigate such planets. In contrast to previous 
studies of habitability in the GL 581 system, we use a consistent atmospheric model to assess surface conditions and habitabil- 
ity. Furthermore, we perform detailed atmospheric simulations for a much larger subset of potential planetary and atmospheric 
scenarios than previously considered. 

Methods. A ID radiative-convective atmosphere model is used to calculate temperature and pressure profiles of model at- 
mospheres, which we assumed to be composed of molecular nitrogen, water, and carbon dioxide. In these calculations, key 
parameters such as surface pressure and CO2 concentration as well as orbital distance and planetary mass are varied. 
Results. Results imply that surface temperatures above freezing could be obtained, independent of the here considered atmo- 
spheric scenarios, at an orbital distance of 0.117 AU. For an orbital distance of 0.146 AU, CO2 concentrations as low as 10 
times the present Earth's value are sufficient to warm the surface above the freezing point of water At 0. 175 AU, only scenarios 
with CO2 concentrations of 5 % and 95 % were found to be habitable. Hence, an additional Super-Earth planet in the GL 581 
system in the previously determined dynamical stability range would be considered a potentially habitable planet. 
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1. Introduction 

The discovery and characterization of extrasolar terrestrial 
planets in the habitable zone (HZ) of their central star is 
one of the most exciting prospects of exoplanetary science. 
Such planets are extremely good candidates for the search 
for extraterrestrial life. The HZ is usually defined as the shell 
around a star where a plane t could retain liquid water on the 
surface dKastin^ etalJll993h . This definition is motivated be- 
cause liquid water seems to be the fundamental requirement 
for life as we know it on Ea r th. Be ing located inside the HZ 
as defined by iKasting et alJ ( Il993l) for an Earth-like planet, 
however, not necessarily implies habitability for a specific 
planetary scenario (see, e.g.. Mars in our own solar system). 
The potential habitability of a planet depends critically on 
atmospheric composition and surface pressure. Still, for a 
planet located well inside this classical HZ, habitability is 
achievable for a much broader range of atmospheric condi- 
tions (greenhouse effect etc.) than for a planet near one of the 
boundaries. 

Among the more than 500 extrasolar planets discovered 
so far, some orbit their c entral star inside or near the 
HZ (e g., Mavor et al. 200 4 iLovis et al.ll2006l iFischer et aP 
I2OO8I lHaghighipouretal.ll2010l) . Most of these planets are 



Neptune- or Jupiter-like gas planets. The planetary system 
Gliese 581 (GL 581), however, contains at least four pl anets 
dBonfils et alJ200 5'.'Udrv et al.'2007'.'Mavor et al.l2009h . one 
of which is a potentially habitable Super-Ea rth, GL 581 d. 
This was shown b y Wordsworth et al. (2010). von Paris et al.l 
(|2010I) . IHu"& Ding C201 1.) and .Kaltenegger et al.. t201 li) who 
presented ID modeling studies of different atmospheric sce- 
narios of GL 581 d. They found habitable surface conditions 
(i.e., surface temperatures above 273 K) with CO2 partial 
pressures as low as 1 bar, depending on CO2 concentration. 
These results imply that the GL 581 planetary system con- 
tains indeed at least one potentially habitable, possibly ter- 
restrial planet. 

Orbita l simulations presented by [Zollinger & Armstrong! 
(l2009l) showed that between the orbits of GL 581 c and 
d (i.e., inside the classical HZ), another Super-Earth planet 
would be dynamically stable. Specifically, they stated a sta- 
bility range for a low-eccentricity planet of not more than 2.6 
Earth masses (m®) ranging from 0. 126 AU to 0. 17 AU. 
Recently, IVogt et al] (l2010l) claimed the detection of two 
more planets in the GL 581 system, one of them (called GL 
581 g in Vogt et al. 2010) with a minimum mass of 3.1 me 
and an orbital distanc e of 0.146 AU, hence inside the stabil- 
ity range calculated bv lZoUinger & Armstrong! ( |2009|) . These 
detections are controversi al and dispute d by further analysis 
of the radial velocity data (lTuomill201 ih . 
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Nevertheless, we use this claimed discovery as a starting 
point to investigate the habitability of planets in the GL 581 
system. Such model calculations aim at supporting the se- 
lection of future targets for detailed observational programs 
of habitable planets wh ich is probably needed in the future 
dHomer & Jonesll2010l) . given the expected number of tar- 
gets. Such potentially habitable planets are expected to be 
discovered in the near future by on-going ground-based pro- 
grams such as MEarth (Nutzman & Charbonneau, 2008) or 



space missions like Kepler (see, e.g., iBorucki et al.ll201 ll for 
an overview of Kepler candidates) and the planned PlaTO 
mission (Catala 2009). First attempts at characterizing the 
atmospheres of transiti ng Super-Earth plane ts have already 
been made (CoRoT-7 b.lGuenther et al .feOl ll a nd GJ 1214 b, 
iBean et alJl20TollDesert et al.ll201 ll lCroll et al.1 l20TTl) . 
For the claimed planet in the HZ of GL 581, dedicated 
modeling studies have been perfo rmed by iPierre humbert 
(|20T1. Heng&Vogt (l20Tlh and Ivon Bloh et all tmuii . 



IPierrehumbert (.201 1.) presented several possible atmospheric 
scenarios (airless planet, pure N2, mixed CO2/H2O atmo- 
spheres) and discussed potential implications for surface 
conditions, without detailed calculations of the atmospheric 
structure for t he mixed CO2/H2O cases. On the other hand, 
iHeng & Vogj (I2OI L) used a general circulation model of 
Earth to simulate the dynamics and circulation on GL 581 
g, however did not investigate surfa ce conditions an d habit- 
ability in detail. The study of v on Bloh et al.l (1201 L) used a 
geodynamic box model to assess planetary habitability, cou- 
pling geophysical and atmospheric processes in a simplified 
approach. 

Previous modeling studies of habitability in the GL 581 sys- 
tem e ither focused on the existing planets GL 581 c and d 
(e.g., ISelsis et all [20071 Ivon Bloh et aP l2007h . used a very 
simple model to simulate atmospheric processes and surface 
conditions (e.g., von Bloh et al. 2011) or investigated only a 
very small subset of potential atmospheric scenarios in terms 
of C O2 level and surface pres sure when varying orbital dis- 
tance dKaltenegger et alj201 lb . We present here model calcu- 
lations for possible terrestrial planets in the GL 581 system 
along the same line of reasoning as in von Paris et al. (2010), 
using a consistent ID atmosphere model. We vary the surface 
pressure and CO2 level over a large range. For these scenar- 
ios, we calculate temperature and pressure profiles in order 
to assess the habitability of up to now hypothetical planets, 
assuming different planetary masses and orbital distances. 
The paper is organized as follows: Sect. |2] states the stellar 
and planetary parameters. The model used is described in 
Sect. [3] A description of the runs is given in Sect. |4] Results 
are described and discussed in Sect. |5] We give our conclu- 
sions in Sect.|6] 



2. Model planets around GL 581 

GL 58 1 is a very quiet M3 star (iBonfils et al.l2005h. The stel - 
lar spectrum of GL 581 is taken from lvon Paris et al. 
It was derived from an UV spectrum measured by the lUE 
(International Ultraviolet Explo rer) satellite and a synthetic 
Nextgen model spectrum (.Hauschildt et al.iri999l) . 
Atmospheric simulations were performed for a subset of 
probable planet scenarios, defined by mass and orbital dis- 



tance. We varied the orbital distance from 0.1 17 to 0.175 AU, 
coveri ng the stability range found bv lZoUinger & ArmstrongI 
(l2009l) . In terms of insolation in the solar system, this trans- 
lates roughly into the present-day insolation at the orbit of 
Mars and the solar flux at Earth about 1.3 billion years ago 
(e.g..l Goughlll98ll ). The planetary mass was varied between 
2.6 and 3.1 m®. For all scenarios, orbits were assumed to be 
circular. 

Th e planetary radius is taken from a mass-radius relationship 
bv lSotin et alJ (l2007h . yielding the surface gravity. Changing 
the planetary mass from 3.1 to 2.6 m^ decreases the gravity 
from 1 6.4 ms~^ to 15. 1 ms~^, i.e by roughly 9 %. 
As in Ivon Paris et al. 1 (l20Toh . the measured Earth surface 
albedo, i.e. the reflectivity of the planetary surface with re- 
spect to incoming stellar r adiation, was taken as the model 
surface albedo (A.„rf^0.13. rRossow & SchiflFeiill99 9') for aU 
scenarios. In doing so, the effect of clouds is explicitly ex- 
cluded from our simulations. Table [T] summarizes the plane- 
tary parameters. Note that scenari o 2 in Table [T]cor responds 
to the claimed planet GL 58 1 g of IVogt et all (l2010h . 



Table 1. Planetary parameters 



Mass [me] Orbital distance [AU] 



Scenario 1 
Scenario 2 
Scenario 3 
Scenario 4 



3.1 
3.1 
3.1 
2.6 



0.117 
0.146 
0.175 
0.146 



3. Atmospheric model 

A cloud-free ID radiative-convective model was used to cal- 
culate the atmospheric structure, i.e. the temperature, water 
and pressure profiles. 

Th e model is originally b ased on the cUmate model described 
bv lKasting et all ( 1198 4a') and' Kasting et al." (1984b). Further 
developments are described by e.g. iKasting (1988) and 
iMischna et af] (j^OO). The model version used here is based 
on the version of ivon Paris et al.l (l2008l) and Ivon Paris et al.l 
(120101) where more details on the model are given. 
The model considers N2, H2O, and CO2 as atmospheric 
species. H2O and CO2 are the two most important green- 
house gases on present Earth, and N2 is present in significant 
amounts in all terrestrial atmospheres of the solar system. 
Temperature profiles from the surface up to a pressure of 6.6 
■ 10"^ bar are calculated by solving the equation of radiative 
transfer and performing convective adjustment, if necessary. 
Convective adjustment means that the lapse rate in the atmo- 
sphere is adjusted to the convective lapse rate instead of using 
the radiative lapse rate if the atmosphere is unstable against 
convection. The convective lapse rate is assumed to be adia- 
batic with contributions of latent heat release by condensing 
water or carbon dioxide. The water pro file is calculated based 
on the relative humidity distribution of iManabe & Wetheraldl 
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(Il967h . Above the cold trap, the water profile is set to an iso- 
profile of the cold trap value. 



4. Atmospheric scenarios 
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We performed a parameter study to investigate the influence 
of surface pressure and CO2 level on the potential habitability 
of terrestrial planets in the GL 581 system, as summarized in 
Table[T] The initial surface pressure (1, 2, 5, 10, 20 bar) and 
CO2 concentration (0.95, 0.05, 3.55 • 10"^ and 355 ppm, re- 
spectively) were varied. N2 was used as the background gas. 
Table |2] summarizes the considered atmospheric scenarios. 



Fig. 1. Temperature-pressure profiles for the low CO2 (355 
ppm CO2) runs of scenario 2. Equilibrium temperature of the 
planet (dashed) and melting temperature of water (dotted) are 
indicated as vertical lines. 



Table 2. Atmospheric scenarios (PAL: Present Atmospheric 
Level) 



The low CO2 runs result in uninhabitable surface conditions 
for all assumed surface pressures (see Fig. [1}, although the 
20 bar run with a surface temperature of 272 K is very close 
to being considered habitable. 



Set 



p [bar] 



Gl (low CO2) 
G2(10PALCO2) 
G3 (medium CO2) 
G4 (high CO2) 



1,2,5,10,20 
1,2,5,10,20 
1,2,5,10,20 
1,2,5,10,20 



CO2 vmr 
3.55 • 10"^ 
3.55 • 10-3 

0.05 

0.95 
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5. Results and discussion 



5.1. Temperature profiles 

The resulting temperature-pressure profiles of scenario 2 in 
Table[T]for the sets G1-G4 are shown in Figs. [T]|4] The equi- 
librium temperature Teq of the planet and the melting tem- 
perature of water (273 K) are indicated by vertical lines. A 
global mean surface temperature of 273 K or higher is gener- 
ally used as the criterion for surface habitability in exoplanet 
science. This criterion is purely based on the phase diagram 
of water where the liquid phase of water needs temperatures 
above 273 K at almost all pressures (the melting line is nearly 
isothermal in the p-T diagram). Note that, of course, on Earth 
life is found in areas with mean annual temperatures far be- 
low the freezing point of water. 



Fig. 2. Temperature-pressure profiles for the 10 PAL CO2 
(3.55 ■ lO"-' CO2) runs of scenario 2. 



When increasing the CO2 content by a factor of 10, calculated 
surface temperatures were ranging between 258 and 320 K 
(Fig.|2]l. For surface pressures of 5 bar or more, surface tem- 
peratures were larger than 273 K, indicating potentially hab- 
itable surface conditions. The rather high surface temperature 
of 320 K for the 20 bar run is the result of the stronger green- 
house effect due to CO2 and a positive water vapor feedback 
(increasing surface temperature leads to increased water va- 
por in the atmosphere, hence more greenhouse effect). 
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Fig. 3. Temperature-pressure profiles for the medium CO2 
(5% CO2) runs of scenario 2. 

The medium CO2 runs (see Fig. O result in habitable con- 
ditions on the surface for pressures of 2 bar and higher 
Calculated surface temperatures are as high as 378 K for the 
20 bar case. 

The high CO2 runs (see Fig. |4]l all show habitable conditions 
on the surface. Calculated surface temperatures range from 
290 to 401 K upon increasing the surface pressure from 1 to 
20 bar. 

In all runs presented here, model atmospheres show convec- 
tive tropospheres and radiative stratospheres, in contrast to 
purely radiative atmospheres encountered in some cases in 
Ivon Paris et al. 1 (12010) . Also, CO2 condensation is absent for 
all runs, even in the high CO2 cases, due to the high temper- 
atures in the middle atmosphere. 
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Fig. 4. Temperature-pressure profiles for the high CO 2 (95% 
CO2) runs of scenario 2. 

5.2. Variation of gravity 

The main effect upon decreasing the gravity g at fixed sur- 
face pressure ps is an increase in column density Dcoi (since 
^^coi ~ This then translates into more greenhouse effect, 
hence the first-order influence of decreasing gravity will be 



an increase in surface temperature. 

However, for the small variations in gravity considered here 
(see Sect. |2]), the resulting increase in surface temperatures 
was only of the order of 0.4 to 2.8 K. Such small variations 
of surface temperatures are not critical in the assessment of 
habitability in the frame of this work and are therefore not 
shown. 

5.3. Variation of orbitai distance 

The immediate effect of increasing the orbital distance is a 
reduction of the stellar flux S received by the planet. Hence, 
the amount of CO2 required to achieve surface temperatures 
above 273 K increases with orbital distance. 
In Table[3] we show the surface temperatures as a function of 
surface pressure and CO2 concentration for the scenarios 1-3 
in Tableffl 

It is clearly seen that for an orbital distance of 0. 1 17 AU even 
the low CO2 runs result in habitable surface conditions over 
the entire range of surface pressures considered. It can also 
be inferred that the same value of surface temperature can 
be achieved for several combinations of surface pressure and 
CO2 concentration. 

As was already demonstrated by the temperature profiles 
shown above, for an orbital distance of 0.146 AU, several 
model scenarios resulted in surface conditions which were 
uninhabitable (e.g., Fig.[T]i. However, in general, habitability 
can be achieved over the entire range of CO2 concentrations 
with accordingly high surface pressures. 



Table 3. Surface temperature in K for variation of orbital dis- 
tance d, surface pressure p in bar and CO2 concentration C. 
Increased grey shading between 260 and 380 K in 30 K steps. 
White indicates cold scenarios. 
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Also towards the outer boundary of the HZ around GL 581, 
at an orbital distance of 0.175 AU, habitable scenarios could 
be found. However, allowed CO2 concentrations for surface 
temperatures above freezing are now limited to the medium 
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and high CO2 case. Model scenarios with less CO2 were 
found to be uninhabitable, independe nt of surface pre s sure. 
These results confirm the findings of Ivon Paris et al. 
that medium CO2 scenarios need to be taken into account 
when assessing the habitability of planets orbiting near the 
outer boundary of the habitable zone. 

6. Conclusions 

We presented ID radiative-convective calculations for a 
subset of potential atmospheric conditions on hypothetical 
Super-Earth planets orbiting GL 58 L We varied parameters 
such as orbital distance, CO2 concentration and surface pres- 
sure. In contrast to previous studies of habitability in the GL 
581 system, we considered a much larger parameter space 
and used a consistent atmospheric model to assess surface 
conditions. 

Our model results imply that habitable surface conditions 
(here Tjurf >273 K) could be obtained for a large part of the 
considered parameter space. For the smallest orbital distance 
of 0.117 AU, habitability was achieved independent of con- 
sidered atmospheric parameter range. For an orbital distance 
of 0.146 AU, depending on surface pressure, CO2 concen- 
trations as low as 10 times the present Earth's value were 
found to be sufficient for surface habitability. However, for 
the largest orbital distance considered (0.175 AU), surface 
conditions were only habitable with CO2 concentrations of 5 
% and more. Hence, our simulations show that an additional 
Super-Earth planet in the G L 581 system in the dynamica l 
stability range calculated by [Zollinger & Armstrongl (l2009l) 
would indeed be considered a potentially habitable planet. 
The model calculations presented here for a subset of the 
possible parameter space (orbital distance, CO2 concentra- 
tion, surface pressure) illustrate how such investigations can 
be helpful to select potentially habitable planets for further, 
more detailed studies from a future larger sample of known 
Super-Earths. 
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